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A system of a two-level atom of an impurity (qubit) inserted into a periodic chain coupled to the 
' continuum is studied with the use of the effective non-Hermitian Hamiltonian. Exact solutions are 

, derived for the quasistationary eigenstates, their complex energies, and transport properties. Due to 

the presence of the qubit, two long-lived states corresponding to the ground and excited states of the 
qubit emerge outside the Bloch energy band. These states remain essentially localized at the qubit 
even in the limit of sufficiently strong coupling between the chain and the environment when the 
super-radiant states are formed. The transmission through the chain is studied as a function of the 
continuum coupling strength and the chain-qubit coupling; the perfect resonance transmission takes 
place through isolated resonances at weak and strong continuum coupling, while the transmission 
is lowered in the intermediate regime. 

. ; PACS numbers: 

^ ' I. INTRODUCTION 
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Open quantum systems are currently in the center of attention of physicists in different subfields. The main driving 
force in this direction is evidently the quest for the nevi^ progress of quantum informatics. Another area with significant 
recent achievements related to openness of a quantum system is nuclear physics, where the understanding of structure 
and reactions of loosely bound nuclei far from stability requires the correct unified treatment of bound states and 
' continuum. Cold atoms in traps and optical lattices can give rise to new effects of coupling, entanglement and transfer 
P3 . of information. Solid-state micro- and nano-devices, including Josephson junctions and spintronics, are probably the 
O most developed arrangements of this type. 

I ^ I , From a general point of view, in all cases we have to deal with a mesoscopic system of interacting constituents that 
serves as a guide for the transmission of a quantum signal. The system can have intrinsic degrees of freedom which 
T— I I can be excited and deexcited by the signal. The coupling to the external world is realized through a certain amount 
^ 1 of channels characterized by the asymptotic quantum numbers of emitted particles or quanta and the final state of 
1^ [ the system. Each channel has an energy threshold where it becomes open and connects to the environment. In the 
: absence of decoherence through external noise or a heat bath, the transmission at given energy is described by the 
' unitary scattering matrix in the space of channels open at this energy. All these features are common for numerous 
• , loosely bound or marginally stable mesoscopic systems determining their main observable properties. 
' A convenient mathematical formalism for description of such systems is given by the effective non-Hermitian 
Hamiltonian; this method based on the Feshbach projection formalism^'^ is formally exact but very fiexible and can 
be adjusted to many specific situations; see the recent review article'^ that covers applications to the continuum shell 
model in nuclear physics and, in less detail, to the quantum signal transmission through simple periodic and disordered 
chains. To stress the breadth of possible problems solved in this approach we can mention recent applications to arrays 
of antennas^ and studies of light harvesting bacteria in the same framework^. 
^ , One bright phenomenon emerging in the situation with a relatively strong continuum coupling in the case when 
the number of open channels is relatively small compared to the number of involved intrinsic states is the so-called 
super-radiance. Being an analog of super-radiance in quantum opticsSi^, this term stands for the formation in the 
system of a collective superposition of the intrinsic states coherently coupled to the same decay channel. The number 
of possible states of this type is equal to the number of open channels. 

The simplest transmission system fully studied with the aid of the effective non-Hermitian Hamiltonian is an open 
periodic chain with hopping between adjacent cells^^— . It was shown^^ that, after appropriate identification of 
parameters, this system is physically equivalent to a realistic sequence of quantum barriers. The transmission was 
studied at weak continuum coupling (a system of narrow individual resonances) and in the limit of strong coupling 
when the resonances overlap and the main role is played by the physics of super-radiance (collectivization through 
continuum). The limit of overlapping resonances is usually described, in the tradition borrowed from nuclear physics, 
in terms of Ericson fluctuations^^; it was shown that this theory should be corrected in a number of aspects, including 
the correct account for the unitarity that is the main source of super-radiance. The behavior of an ideal periodic 
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system was juxtaposed to the case of disorder, where the super-radiance still survives, and it is possible to establish 
the connection to the Anderson model, degree of chaos inside the system and universal conductance fluctuations. The 
consideration was also extended to grids with two- and three-dimensional geometry; an especially interesting case is 
presented by the star graph^ where a number of open channels intersect at a common central point. At this point 
an analog of a bound state (evanescent wave) with a long lifetime exists that can serve for accumulation of quantum 
information. 

Thanks to the development of nanotechnology, the preparation of low-dimensional assemblies of nanoparticles 
becomes a routine experimental task^^. This renews the attention to the study of quantum properties of relatively 
simple low-dimensional mesoscopic systems which reveal a rich physical behavior. In this context, the exactly solvable 
models as tight-binding one- or two-dimensional chains coupled to adatoms or impurities are of special interesti^^— . 

In what follows we consider a similar system with a two-level atom of an impurity (qubit) inserted into an open 
periodic chain. The model is simple enough to allow for the exact solution; at the same time the model turns out to be 
rich enough to demonstrate interesting physics of the signal transmission. An analytical consideration is supplemented 
by the detailed numerical study for the chain with eleven sites. 

The paper is organized as follows. In Sec. 2 we consider a closed chain of 2N identical cells with the nearest 
neighbor hopping interaction. Two arms of the chain are connected through the central cell occupied by a two-level 
atom (qubit). We find the energy spectrum of the system that consists of a normal band of delocalized Bloch standing 
waves and two additional states outside the band corresponding to the excited and the ground states of the qubit. 
We study the evolution of the energy spectrum as a function of the chain-qubit coupling strength. 

Sec. 3 describes the energy spectrum of the open system coupled to the continuum through its edge states. The 
former stationary states acquire decay widths that change as a function of the coupling constants. The decay widths 
of the qubit states are small compared to those of the Bloch waves. The qubit states remain essentially localized at 
the qubit even when the coupling of the chain to the continuum becomes sufficiently strong to allow for the formation 
of super-radiant states. 

In Sec. 4 we study the transmission through the chain as a function of the coupling parameters. In both limits of 
the edge continuum coupling being weak and strong compared to the coupling between the qubit levels, the resonances 
are well isolated all having narrow decay widths and perfect transmission at resonance energies. When both couplings 
are of comparable strength, the resonances are overlapped and the transmission is below the perfect level. The results 
are summarized in the Conclusion. 



II. CLOSED CHAIN 

We consider a linear chain of 2A'' identical cells numbered asn — —N, —{N~l), —1 and n = 1, ...,N — 1, N, while 
the central cell, n = 0, is occupied by a qubit, a two-level atom with states |0) and |e), excitation energy eg — eo = A, 
and matrix element A of the qubit excitation. The energy eg is the level position in all cells; for simplicity we put the 
lower level of the isolated qubit at the same position. Introducing the hopping (tunneling) matrix element v between 
the neighboring cells, we come to the Hamiltonian of the closed chain: 

i/„„ = eo, Hn,n+1 = Hn+i^n = V, n = -iV, 0, iV; (1) 

Hee — A, Hoe — iJoe = A, (2) 

where all matrix elements can be considered as real. 

Before introducing the coupling to the outside world, we briefly characterize the solution for the closed chain. A 
general stationary state \E) with energy E can be presented as a superposition 

N 

\E)^ J2 c„{E)\n)+b{E)\e). (3) 

n=-N 

The boundary conditions for the chain closed at the edges are c_jv-i ~ cn+i = 0. The coefficients of the superposition 
satisfy the obvious equations: 

{E - eo)cn - v(cn-i + c„+i) ^ 0, n^O, (4) 



{E - eo)co - v{c^i + c+i) = Xb, 



(5) 
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{E - A)b = Acq. (6) 
Eqs. (O and © can be treated as boundary conditions for the two chains (left and right) implying that 

Following the same procedure as for a single chains', see also Appendix in Ref J^, we find that the solutions on both 
sides of the qubit have the form 



A'^l + B'C, n>0. ^ ' 



The amplitudes ^± are given by 



E-eo±^/(E-eoy-4v^ , (9) 



and ^+^~ = 1. The edge conditions determine 



B = -A^^-j , B' = -^'(^i±j . (10) 

The central point determines two classes of solutions. The antisymmetric states, ci + c_i — 0, have cq = 0, so that 
these wave functions are decoupled from the excited qubit that lives on the excited level i? = A (in this situation the 
qubit can be excited or deexcited only by an additional external coupling). Because of cq = and interaction only 
between the neighboring cells, the two sides of the chain are also decoupled from each other so that such modes do 
not take part in the transport through the whole chain. The spectrum Eq of the antisymmetric states is determined 
by the roots of = i that can be parameterized by the even-number quantized quasimomentum q, or by the 

phase (fiq, 

Eq = 2vcosipq, (fq^ 2 TV + 2 ' ^^^^"^^ (11) 
and the local amplitudes ([8|) of the wave functions are 



The energies (fTT|) are inside the Bloch band (— 2w, +2v). 

The symmetric solutions, ci = c_i, involve the qubit into dynamics and the corresponding roots are determined by 



.2^+2 ^ + {E - A)^E - e,^^ - Av^ 
^+ X'-{E-A)^{E-eoY^Av^' 

To derive this relation, it is convenient to use eq. (|4]) and connect the amplitude ci, or c_i, to the central amplitude 
Co, 

With the qubit disconnected, A 0, we have the missing in eq. (jlip g-odd part of the band spectrum with the 
amplitudes ()12p where sm{nipq) is changed to cos{rnpq). When the qubit is connected to the chain, the upper level 
is outside the band, E > 2v; the lower level also leaves the band at a finite value of A when E < —2v. The full 
spectrum of twelve energy levels for = 5 is shown as a function of A in Fig. 1. The upper state starts at i? = A and 
grows approximately linearly with A. At small A, this is the state with the excited qubit and only weak admixtures 
of intrinsic sites. 

With increasing A, the symmetric and antisymmetric states inside the band become degenerate, while the excited 
(upper) wave function is spread almost equally over two states of the qubit. The orthogonal combination of the 
excited and ground state of the qubit gives the lower state (at sufficiently large values of A). In this limit the picture 
effectively is of Rabi oscillations between the qubit levels with only a small probability of hopping along the chain. 
The chain then is almost decoupled from the qubit. This evolution of the two wave functions is illustrated by Fig. 2. 
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FIG. 1: (Color online). Energy levels for a system of the closed chain of eleven sites and the excited qubit state in the middle 
as a function of the qubit coupling strength A; the hopping amplitude is set to n = 1, and the excitation energy of the qubit 
A = 2.5. 




FIG. 2: (Color online). Squared components of the lower (left column) and upper (right column) of the wave functions as a 
function of the qubit excitation strength A, for a chain N = 5 and i; = 1. 
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III. OPEN CHAIN 

Now we assume that the edge states are coupled to the outside world by ideal leads and the situation becomes 
identical with what can be described by the effective non-Hermitian Hamiltonian, 

n = H-'-W, (15) 

acting only in the intrinsic space of the closed system. Here the anti-Hermitian part W is factorized- having the 
matrix elements W12 between the intrinsic states |1) and |2) in the form of the product of partial amplitudes Al and 
A2 which generate the interaction between intrinsic states through an open decay channel c, 

Wi2^ AlA^. (16) 

c; open 

It is important that here the superscript c runs only over open channels {on-shcU interaction). In general the amplitudes 
Ai are energy-dependent vanishing at the threshold energy for a given channel, and their near-threshold behavior 
determines the non-exponential decay curve in the long-time limit. In many cases this dependence can be ignored, 
which we assume in this work as well. 

In the case of a linear chain we allow only left and right decays through the edge states (it would be also interesting 
to study the decoherence through the coupling with many "random" weak channels connected to the intrinsic sites). 
The factorized nature of the operator W (dictated essentially by requirements of unitarity of the scattering matrix 
in the channel space, see for example^^), shows that this operator has only few non-zero eigenvalues their number 
being equal to the number of open channels. The corresponding eigenstates are obviously just the edge states directly 
coupled to the continuum. Therefore it is sufficient in our scheme to introduce two complex energies, and for the 
states at the edges, where eL,R = eo — {i/2)"fL,R- Such an open system without a qubit was analyzed in Refs.— >^'i^'i^. 
In this approach the states directly coupled to the continuum play the role of doorways^S., and the remaining states 
can get their widths (finite lifetimes) only through their coupling to the doorways. 

The typical situation for the chain without a qubit is evolving as a function of parameters ^l,r compared to the 
level spacing D in the closed system. At weak coupling, every intrinsic state becomes a resonance with a small decay 
width determined by the overlap of the Bloch state with the edges; the final width distribution has, for 7^ = 7^^, a 
maximum in the center of the band. In the limit of strong continuum coupling, we have a super-radiant situation when 
the central state accumulates almost the entire width while the remaining states become very long-lived (trapped). 
In the case of 7l ^ jr there occur two super-radiant transitions^^-- with the maximum of the signal transmission in 
between. In the site representation, the super-radiant states with energies in the center of the band are located at the 
edges. This picture survives also the possible presence of disorder in the intrinsic wells^". 

The results for the chain with the qubit are illustrated by the series of graphs, where the chain consists of 5-1-1-1-5=11 
cells with the qubit in the middle, altogether 12 intrinsic states. Here we diagonalize the effective Hamiltonian in the 
doorway representation: the continuum coupling occurs only at the edges which serve as doorways, and the matrix 
elements of the anti-Hermitian part of the effective Hamiltonian which couples the states \q) with the outside 
world, are given in the band representation by 

Wqq' = ^LC^N{q)c^N[q') + -lRCN{q)cN{q')- (17) 

Fig. 3 shows the evolution (as a function of A) of the resonance energies for the case of weak continuum coupling, 
7l = 7fl = 0.1 (the scale is fixed by the band width, v = 1). At small A, we have the parabolic distribution of 
widths with the maximum at the center of the band, as known from previous studies^iiS, and the decoupled excited 
qubit state above the band with zero width. As A increases, the symmetric and antisymmetric states merge becoming 
effectively decoupled from the qubit. The excited qubit state and emerging outside the band the ground qubit state 
are still almost stationary (Rabi regime). They are moving along the real energy axis being repelled by the band. 

The situation changes when we come to the strong continuum coupling. Fig. 4. Here we again follow the evolution 
as a function of A but at 7 = 20. At small A = 0.1, the two coinciding super-radiant states are formed in the center of 
the band (here we keep — 1r), and the remaining ten states are trapped, effectively returning to the non-overlap 
regime (the left lower plot shows, at a much smaller width scale, the parabolic width distribution). With increase of 
A, the super-radiant states survive, while the qubit states again are repelled by the band along the real energy axis 
having still very small widths. 

It is instructive to take a look of the width evolution as a function of the continuum coupling strength 7. The 
typical process is presented by Fig. 5, where the trajectories of all twelve complex poles in the lower half of the 
complex plane are shown. In the limit of very weak continuum coupling, part (a) of this figure, the width distribution 
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FIG. 3: (Color online). Resonance complex energies (eigenvalues of the effective Hamiltonian for weak continuum coupling, 
7 = 0.1) evolve as a function of the qubit excitation strength A. Two qubit states are effectively decoupled from the chain and 
move along the real energy axis. 
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FIG. 4: (Color online). Resonance complex energies (eigenvalues of the effective Hamiltonian for strong continuum coupling, 
7 — 20) show segregation of two super-radiant states in the middle of the band from ten trapped states which include the 
strongly localized states of the qubit. 
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FIG. 5: (Color online). Complex-plane trajectories of eigenstates of the effective Hamiltonian; the parameter values are 
V = 1, \ — 0.8, A = 2.5. Panel (a) shows the behavior of 12 states in the chain while panel (b) singles out the state genetically 
connected to the excited qubit state located outside the band. 

is parabolic, being proportional, as a function of real energy, to the group velocity of the band states. It is transformed 
with increase of 7. All widths, except for the super-radiant states in the center of the energy band, turn back after 
reaching their maximum values. Their corresponding trajectories are almost symmetric with respect to their maxima 
which is typical for the phenomenon of super-radiance in a space of fixed dimension. Indeed, after the segregation 
of the super-radiant state(s), the remaining trapped states are essentially in the same situation as they were in the 
beginning of the process; this symmetry for the finite dimension of intrinsic space is a characteristic featurei^ that 
appears also in the statistical distribution of neutron widths for thermal-energy neutron resonances^. Part (b) of 
the figure selects, on a detailed energy scale, the complex-plane evolution of the excited qubit state that becomes 
extremely long-lived. The details of interference between neighboring resonances were also discussed repeatedly in 
the context of the electron conduction in nano-scale systems, see for example-^. 

It was noticed long agc^^ that, in the description of an open system with the aid of the effective Hamiltonian, the 
Hermitian and non-Hermitian parts of the interaction act in the opposite way. The real (Hermitian) perturbation 
repels the levels but, through the mixing mechanism, attracts the width of unstable states. Contrary to that, the 
imaginary (non-Hermitian) interaction through the continuum repels the widths (the road to superradiance and 
trapping) but attracts real energies of resonances. This attraction is seen in Fig. 6 for a larger value of A: as the 
continuum coupling 7 increases along the road to super-radiance, the real energies of the poles move to the middle of 
the band, where the super-radiant states are located. 

Finally, the limit of very strong coupling between the qubit levels is shown in Fig. 7 where A = 8. As antisymmetric 
and symmetric roots merge with increasing continuum coupling 7, we see, panel (a), the evolution of pairs of states 
(ten of them inside the band including the super-radiant ones). The time arrow of this evolution is again in the 
direction of attraction for the real energies of resonances. The qubit states, the ground state, (b), and the excited 
state, (c), are essentially decoupled. They do not shift and only gradually increase their (still small) decay widths as 
7 increases. 
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FIG. 6: (Color online). Complex-plane trajectories of eigenstates of the effective Hamiltonian; the parameters values are 
V = 1, X — 2.1, A = 2.5. Ten states, including the super-radiant at the edges of the chain, are moved closer inside the band, 
panel (a), two states mainly localized at the qubit have energies outside the band and essentially interact only with each other 
having a large lifetime with respect to tunneling through the chain, panels (b) and (c). 
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FIG. 7: (Color online). Complex-plane trajectories of eigenstates of the effective Hamiltonian; the parameters values are 
V = 1, A = 8, A = 2.5. Ten states, including the super-radiant at the edges of the chain, are moved closer inside the band, 
panel (a). The qubit states, the ground state, (b), and the excited state, (c), are essentially decoupled. 
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IV. TRANSMISSION THROUGH THE CHAIN 



Here we briefly consider the transmission of an external signal through the chain with the inserted qubit, similarly 
to the consideration made earlier for the uniform chain or multi-dimensional lattices^^ and for the star graphic. The 
internal propagation of the signal of given energy E in the open system is described by the propagator 

CiE) . (18) 

where Ti is the effective Hamiltonian (15) with the imaginary part (|16p that describes multiple excursions of the signal 
into continuum and back. Our simple geometry has two open channels, left and right (labeled L and R). The full 
amplitude Z'"^ of the process between channels a and b starts with the entrance amplitude A^^ that populates the 
intrinsic state \n) and ends with the exit amplitude A^^ from the intrinsic state |m); all paths b ^ a interfere: 

Z'"^{E)^Y.'^'nGmn{E)An- (19) 

mn 

It is easy to show that the corresponding scattering matrix, S^'^ — 6^°" —iZ^'\ is unitary^.. The transmission coefficient 
is given by 

T^°-{E) ^\Z^''{E)\^ . (20) 

It might be convenient to perform the transformation to the (biorthogonal) basis \r) of eigenfunctions of the effective 
Hamiltonian. The complex energies £r — Er — {i/2)Tj. correspond to the poles of the scattering matrix, while the 
process amplitude (|19l) still has factorized residues transformed to the eigenbasis. 

We can note parenthetically that this description can be treated as a simple superposition of interfering resonances 
only approximately, namely if the energy dependence of continuum amplitudes A'^ is neglected as it is done in our 
consideration. In this approximation the time decay curve of a single isolated resonance \r) would be pure exponential 
with the width F^. 

For the calculation of transmission we adopt the approach of Ref.^^. For the open chain of Sec. 3, the transmission 
is determined by the edge couplings which we again assume here to be equal, Jl ^ 7r = 7- Similarly to Ref.^^, the 
transmission coefficient (l20l) can be written as 



rpRL _ rpLR. 



hy)iE - A) 



(22) 



Below we show the transmission results for the chain of iV = 5 (twelve intrinsic states) and various combinations of 
the parameters. The resulting picture is determined by the counterplay of the trend to super-radiation and decoupling 
of the qubit. 

Starting with the weak qubit excitation amplitude. Fig. 8 for A = 0.1, we follow the evolution of the transmission 
as a function of the continuum coupling 7. At small 7 = 0.1, panel (a), we see twelve isolated resonances all having 
narrow decay widths; one of them is outside the energy band as we discussed earlier. When 7 is growing, panels (b) 
and (c), the resonances start overlapping. When two super-radiant states merge at 7 = 2.4, panel (c), they disappear 
from the transmission spectrum, so that panel (d) for 7 = 4 shows ten separated resonances. At each resonance the 
transmission is perfect, T — 1. 

Next four panels. Fig. 9, correspond to the intermediate value A = 2. Again the case of weak continuum coupling, 
panel (a), 7 — 0.1, reveals twelve resonances; now the two states associated with the qubit are outside the energy 
band, while all resonances still show perfect transmission. At 7 = A = 2, panel (b), the competition between the 
transmission through the chain and Rabi dynamics of the qubit leads to an almost random pattern of overlapping 
resonances with transmission below perfect, similar to Ericson fluctuationsiS or universal conductance fluctuations^—. 
The relation between those well known pictures and necessary changes due to the effects of unitarity in exact theory 
were discussed in Ref.— . At large 7 ^ A, panels (c) and (d), the continuum coupling prevails leading to the narrow 
resonances coming from trapped intrinsic states. 

Finally, the large value of the qubit excitation strength, A = 5, changes the transmission picture at not very large 
7, Fig. 10. At small continuum coupling, 7 <C A, panel (a), the two states associated with the qubit produce two 
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FIG. 8: (Color online). Evolution of the transmission for small A = 0.1 as a function of the continuum coupling parameter 
7l = 7fl = 7, from 12 isolated resonances (one outside the energy band) through overlap and super-radiance to ten resonances 
corresponding to trapped states. As earlier, A — 2.5, v — 1. 




FIG. 9: (Color online). The same as Fig. 8, with A = 2. 
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FIG. 10: (Color online). The same as Fig. 8, with A = 5. 



resonances outside the band, with transmission equal to 1. At 7 = 2, panel (b), we observe some kind of an intrinsic 
resonance between propagation and internal oscillations which, along with the emergence of super-radiance, almost 
kills the transmission at other energies within the band. With further growth of 7, panels (c) and (d), the perfect 
transmission through trapped states, including those associated with the qubit, is gradually restored. This abundance 
of possible regimes opens the way to various applications. 



V. CONCLUSION 



Using the formalism of the effective non-Hermitian Hamiltonian, we studied a model of quantum signal transmission 
through a linear periodic chain with a qubit placed at the center. For the closed chain, the intrinsic eigenstates form two 
classes, antisymmetric standing waves with the excluded center site and two independent subchains, and symmetric 
that includes the qubit dynamics with two additional eigenstates. When the system is coupled to the environment, 
we have found the spectrum of quasistationary states characterized by complex energies and finite lifetimes. Two of 
those states are genetically related to the ground and excited state of the qubit with real energy outside the Bloch 
band. 

The most interesting feature is the stability of states related to the qubit, small decay widths and correspondingly 
long lifetimes. Further, these states are only weakly perturbed when the coupling of the chain to the continuum 
becomes strong. In the limit of strong coupling, the qubit states and remaining trapped Bloch states are practically 
shielded from the external world by the two super-radiant states localized at the edges of the chain. The stability 
of the qubit states and the possibility to switch on and off the access to them suggests that the simple configuration 
considered above may serve as a building block of a quantum computer. 

We discussed also the transmission through the chain as a function of the coupling strength to the continuum. In 
the limits of weak and strong continuum coupling (as compared to the excitation amplitude of the qubit) the chain 
reveals well separated narrow resonances with perfect transmission at corresponding energy. In the intermediate 
regime, when the continuum coupling and the excitation strength of the qubit are comparable, the resonances overlap 
with the transmission below the perfect level. 

There are many possibilities to enrich this prototypical model. For a realistic situation, for example a chain of 
quantum dots, one should carefully determine the lifetimes of the qubit states. The geometry of the system can be 
made more complicated in various ways including the transition to more-dimensional schemes. More qubits and more 
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branches can be added approaching a comphcated network. It would be also interesting to extend the ideology of an 
open quantum system to the study of coherent photon transport in circuit quantum electrodynamics^'— ■2?^. 
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